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Abstract 

The branching ratios for the hadronic B-meson weak decays B — > J/ipK and 
B — > Dtt are used to extract the size of the "non-factorizable" terms in the decay 
amplitudes. It is pointed out that the solutions are not uniquely determined. In the 
B — > J/i/jK case, a 2-fold ambiguity can be removed by analyzing the contribution 
of this decay to B — > Kl + l~ . In the B —* Dtt case, a 4- fold ambiguity can only 
be removed if the "non-factorizable" terms are assumed to be a small correction 
to the vacuum insertion result. 
PACS:13.25.Hw,13.25.-k,12.15.Ji 



1 Introduction 



An increasing sample of 5-mesons has been gathered from different exper- 
iments in recent times, and will tend to increase sharply in the near future 
with the advent of 5-factories and, possibly, experiments in hadron colliders 
targeted at 5-physics. The major concern in this paper is to clarify how 
such wealth of data can be used to study some of the aspects that remain 
unclear in the hadronic weak decays of the B and the other flavored mesons. 
The focus shall be on the 2-body decays that proceed through the tree level, 
Cabibbo favored, quark transitions b — > ccs and b — ► cud. The correspond- 
ing effective weak Hamiltonian, once QCD corrections have been included, 
is the sum of two 4-quark operators, that only differ in the color indices of 
their quark fields and the strength of their Wilson coefficients. In the calcu- 
lation of the decay amplitudes, one is faced with the task of evaluating the 
matrix elements, between the initial and final hadronic states, of those two 
4-quark operators. The vacuum insertion (factorization) approximation [|TJ 
reduces this problem to that of determining the matrix elements of bilinear 
quark operators; such matrix elements can then be measured from leptonic 
and semi-leptonic decays, or calculated in some model for the mesons. Un- 
fortunately, this is a poor approximation for one of the 4-quark operators: 
due to a mismatch in the color indices, its matrix element has significant 
"non-factorizable" terms [0, that do not appear in the vacuum insertion ap- 
proximation. That these terms are important can be seen, for example, in 
the strong disagreement between the factorization predictions and the ob- 
served rates for the color suppressed B or D decays (see, for example, ref. 
[§]. In the color suppressed decays, the effect of the "non-factorizable" terms 
is enhanced by an accidental cancellation of the terms that are non-zero in 
the vacuum insertion approximation). 

The standard procedure in dealing with this discrepancy has been to 
preserve the vacuum insertion result for the hadronic matrix elements, but 
to replace the Wilson coefficients that multiply them by two free parameters. 
These parameters are then determined from a fit to the observed values of 
the branching ratios. For the case of the D- and 5-meson decays, the two 
parameters (for each case) fit the available data quite well. For the D-mesons, 
the values of the parameters correspond to dropping the contribution of the 
color mismatched 4-quark operator altogether , which is a procedure that 
finds some theoretical justification in l/N c expansion arguments 0. Quite 
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surprisingly, the recent data on 5-meson decays has shown that, in this case, 
the values of the parameters do not obey the same rule Q: neglecting the 
contribution of the color mismatched operator cannot be used as a systematic 
procedure to obtain the value of the free parameters, as it was the case for 
the D decays. 

In view of the failure of the factorization approximation when applied to 
the color mismatched operator, and the failure of the standard procedure 
of dropping the contribution of that operator altogether, when applied to 
the 5-meson decays, the "non-factorizable" terms have to be dealt with. 
Theoretical estimates, based on QCD sum rules, have been presented in 
refs. H and 0, for some decays. Here, I address the more basic question 
of extracting the size of the "non-factorizable" terms from the data. This 
program was first advocated by Deshpande, Gronau and Sutherland 0, and 
it has been recently applied to both D- and 5-meson decays by Cheng M. 
In this work, I concentrate on the case of the B decays; where, except for 
the matrix elements of the operators with a color mismatch, the vacuum 
insertion approximation can be assumed to work well || (in particular, the 
effects of inelastic final state scatterings will be neglected for the 5-decays 
considered in here). Special attention is paid to the ambiguities in the values 
of the "non-factorizable" terms that are extracted from the data. The results 
of Cheng M are recovered among other possible solutions. Finally, ways of 
lifting these ambiguities are discussed. 

2 The "non-factorizable" terms in B —> J/ipK 
and B —> Dn 

The tree level, Cabibbo favored, hadronic weak decays of the 5-mesons cor- 
respond to the quark transitions b — > ccs and b — * cud. The decay amplitudes 
are derived from the effective weak Hamiltonian 

Heff = ^ WdYZfrMOZ + C 2 (n)O c 2 ) 

+V A V: d {C x { l i)Ot + C 2 { l i)Ot)], (1) 

where 

0\ = c a ^{l - 7 5 )6 Q 3/37^(1 - 75)c/3, 
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0% = ^(1-75)^5^(1-75)^, (2) 

and the operators 0\ 2 are obtained from Of 2 > replacing s and c by d and 
w, respectively. The Wilson coefficients Ci(/x) and C 2 (//) contain the short 
distance QCD corrections. In the leading logarithm approximation, they are 
Ci >2 = (C+ ± C_)/2, with 

^M-f^r)*""' (3) 

(7_ = — 27+ = 2; rij is the number of active flavors). For = 200 MeV 
[fLlfl , and at the scale fi = 5.0 GeV, this gives 

Ci = 1.117 C 2 = -0.266. (4) 

For the exclusive decays B~ — > J/ipK~ or B® — > J/ipK°, the amplitude 

is 

= -^K 6 V^(d < J/^K\0{\B > +C 2 < J/^K\O c 2 \B >). (5) 

The hadronic matrix element of 2 is calculated in the vacuum insertion 
(factorization) approximation: 

< J/^K\O c 2 \B >=< 7/^15^(1-75)^10 >< K\s al »{l- lh )b a \B > . (6) 

Whereas, for O^, an additional "non-factorizable" term is included that re- 
flects the color mismatch in the c-c quark fields: 

< J/ipK\<D{\B >= 

= — < J/^c^l - 75)c/?|0 >< K\s a ^{l - ^)b a \B > 

+ < J/^K\0{\B > non _ /act . . (7) 

The factor 1/N C corresponds to the projection of the two color mismatched 
quark fields into a color singlet; it is included so that the "non-factorizable" 



3 



term (that is defined by eq. [7]) vanishes in the vacuum insertion approxima- 
tion. The decay amplitude can then be written as 

A{B - J/4K) = -^V cb V; s Ma, (8) 

with 

a = C 2 + C 1 (^ f +X), (9) 

M =< J/^\c plll {l - 75)^10 >< K\s a ^(l - 75 )6 a |S > (10) 

and 

X EE < J/^K > non _ /act . . (11) 

In the BSW model @, M = 5.84 GeV 3 x /^/(395 MeV). With \V cb \ = 



0.038^1.63 psec/r B |T3|, and = 395 MeV (which corresponds to T(J/ifj 
e + e~) = (5.26 ± 0.37) keV fLlp, the branching ratio is 



B(B -> J/^K) = 1.90|a| 2 %. (12) 



From the average of the experimental results [11 



B(B~ — > J/ipK~) = (0.102 ±0.014)% 
B(W ^ J/ipW) = (0.075 ±0.021)%, (13) 

it follows that \a\ ~ 0.22, and so 

X ~ -0.29 or 0.10 (14) 

(the numerical differences with respect to the analogous results in ref. || 
correspond to the updated values of the parameters that were used in here). 
The 2-fold ambiguity in the value of X, which corresponds to the unknown 
sign of a, cannot be resolved by the branching ratio of B — > J/ipK alone. 

Proceeding in a similar way for the exclusive processes B% — > D + tt~, 
B® — ► D°7T , and B~ — > D°7r~, the weak decay amplitudes are 



,00 

V2' C °'" d v / 2 J 



^4 + ~ = ^K*K^M iai3 



A 00 = —Lv cb V: d —M 2 a 2 
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and 



lit— (•'! ,- Tr * r, , M 2 a 2l 



f- = ^ TO f iai [l + ^^ (15) 

respectively (non-spectator contributions are very small and they have been 
neglected). The hadronic matrix elements Mi and M 2 are 

Mi =< £> + |c a 7^(l - 7 5 )M#° ><; ^1^(1 _ 75 ) U/3 )|0 > (16) 

and 

M 2 = V2< ir°\d a ^(l - 75)&a|^ D |c^(l - 75 H|0 > . (17) 

In the BSW model 0, M x = 1.85 GeV 3 and M 2 = 2.28 GeV 3 (for f D = 220 
MeV). As in the previous case, the parameters 

a 1 = C 1 + C 2 (^ + X 1 ) (18) 



and 



include the terms 



a 2 = C 2 + C 1 (±r+X 2 ) (19) 



1 

/ \~r 1 /' 1 " 1 / j , , 

c/ ^non-fact. 



Mi 1 2 



X 2 = ^ < D°7r°|<^|SS > non -fact, (20) 
where the "non-factorizable" matrix elements are defined by 



< D + ti-\0^\B° > = —M l +<D + Ti-\O%\B d > 

1 



2 I d -^non-fact. 



V2<D\ \O^\B° d > = — M 2 + V2<D 7r°\Ol\B° d > non _ fact , . (21) 

In order to determine the values of the parameters ai and a 2 , the magni- 
tudes of the decay amplitudes are extracted from the experimental value of 
the corresponding branching ratios JET], 

B(B% -> D+tT) = (0.30 ±0.04)% 
B(W d ^D°n ) < 0.048% (90%(7.L.) 
B(B-^D\-) = (0.53 ±0.05)%, (22) 
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and compared to the predictions in eq. ISJ. The later are the amplitudes in 
the absence of final state interaction effects, and so the comparison must be 
done with care (to be sure, I use the notation A + ~ , A 00 and A°~ for the full 
amplitudes). 

As for the B — > J / ipK case, it is assumed that the effects of inelastic final 
state interaction scatterings are negligible. However, the B — > Drr decays 
involve two isospin channels, and so an elastic final state interaction phase, 
5, can appear between the two isospin amplitudes A 3 / 2 and Ai/ 2 . In general, 

A3/2 = 1^4-3/2 1 

A 1/2 = \A 1/2 \e iS ', (23) 

where 5' = 5 or 5 + tt (according to the relative sign of the two amplitudes 
in the absence of final state interactions). The full amplitudes A + ~ , ^4 00 and 
A°~ are related to the isospin amplitudes A 3 / 2 and A 1 / 2 , in the following 
way: 

A+ ~ = TS^ + S- 4 " 2 ' 

A*® ^ A 1 A 

A = -^V--^, 

A ' = V3A 3/2 . (24) 
This allows to determine the magnitudes \Ai/ 2 \ and |v4 3 / 2 |, as well as cos 5', 



from the experimental results in eq. 22. In particular, it follows that 



cos5'> 0.77 (25) 

The lack of a more precise value for 5' is due to the fact that only an upper- 
limit exists for B(B% — > D n°). For now, I will take 5' = (i.e. the final 
state interaction phase is either 5 = or 71). Then, 

\AiM = ^(l^l-^I^D- ( 26 ) 
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Since the decay amplitudes in the absence of final state interactions are 

A°° = ^\A 3/2 \^l=\A 1/2 \, 



A°~ = Vs\A 3/2 \ (27) 



(the two signs correspond to 8' = 8 or 8 + ir), it follows that 

A°~ \A 
A+- ~ \A 

for 5 = 0: and 



A°- \A°-\ 

\A + -\ = \A+-\ 4zr = ]4zzk (28) 



\A + -\ = \ 2 -\A°-\-\A+-\\ ^= 2l J_? (29) 



\A°- 




2 
3 







for 8 = 7r. The predictions of eq. [15], in terms of the parameters a± and a 2 , 
are replaced on the LHS of these equations; whereas the experimental input 
from the branching ratios is used for the RHS. If the experimental result of 



eq. 25 is interpreted as showing a negligible phase shift from the final state 



interaction effects, eq. |8| gives (with \V c b\ as before) 



lad ~ 1.07 1 + 1.23— ~ 1.33; (30) 

CLi 

and the size of the "non-factorizable" terms is then 

X x ~ -0.17 X 2 ~0.16 (31) 

or 

X x ~ 7.90 X 2 ~ -0.35 (32) 

(the positive value for X 2 corresponds to the result in ref. ||). Alternatively, 
the data can be interpreted as showing a maximal phase shift. Then eq. |2~9| 
gives 

|ai| ~ 0.12 1 + 1.23— ~ -11.63; (33) 

ai 
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and so 

X 1 ~ 1.70 X 2 ~ -0.61 (34) 

or 

Xx ~ 2.16 X 2 ~ 0.52. (35) 

The 4-fold ambiguity corresponds to the fact that the final state interac- 
tion phase can only be determined modulo n, and the sign of a% cannot be 



determined from the branching ratios in eq. |22. 

At this point, a word should be said about the uncertainties in the re- 
sults that were presented. The derivation of the parameters \a\, \ai\ and \a 2 \ 
suffers from the experimental errors in the branching ratios (in particular, 
B(B^ — > D°tt°) is still missing), and in \V c b\- These will improve with more 
accumulated data; which will also allow to derive the hadronic matrix ele- 
ments of the bilinear quark operators from the semileptonic branching ratios 
(see the tests of factorization in ref. for example). At present, the use 



of the BSW model [12| entails an uncertainty that is hard to quantify. The 
derivation of the terms X, X\ and X 2 suffers from the additional uncertainty 
on the Wilson coefficients: in particular, the value chosen for the scale fi is 



important |L4] . The positive solution for X and the negative solution for X\ 
are the most sensitive as they involve somewhat delicate cancellations; for fi 
in the range 4.5 to 5.5 GeV, they oscillate by as much as 25% (whereas the 
other results vary by not more than 10%). I have taken /i = 5.0 GeV, which 
is approximately the constituent 6-quark mass. 

It must be stressed that the vacuum insertion result has been assumed 
to work well for the operators with the correct color assignments (see for 
example eq. ^), and the inelastic final state scatterings have been deemed 
negligible. This is necessary in order to be able to derive the size of the 
"non-factorizable" terms from the data. In ref. || , it is argued, on the basis 
of color transparency, that these assumptions are expected to hold for the 
low-multiplicity B decays. The argument is less reliable for the case of the 
D decays, and so they have not been considered in here. 



3 Resolving the ambiguities 
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3.1 "Non-factorizable" terms in B — > J/ipK 

The sign of the parameter a, that appears in the B — > J '/ipK amplitude of 
eq. II can be determined from the interference between the short distance 
contribution to B — > Kl + l~ , and the long distance contribution due to B — > 
KJ/ip —>■ Kl + l~ . The short distance amplitude is derived from the effective 
weak Hamiltonian in eq. [I] (the operators 0\ 2 contribute at the 1-loop level), 
with the additional electroweak terms: 



Hi 



eff 



%v th vz £ cMd 

i=7,8,9 



(36) 



where the operators 



o 7 
o 9 



a 



;m b s a a" v (l + 'y 5 )b a F t 



[IV 



— s„7 M (l - 75)W7/J 

271 



a 



— s a 7 M (l - 7 5 )W7m75^ 



contribute to B — > Kl + l at tree level. For m t 



(37) 

175 GeV, the Wilson 



coefficients in eq. |36], in the leading logarithm approximation, are [ff7| C-j 
0.326, C$ = —3.752 and Cq = 4.581. The decay amplitude is 



G F a 1 

+ < 7^(7^(1 + -f 5 )b\B > C 7 m b ^ua^vj\ 
= Pb — Vk)- The factor 



(38) 



a 



Seff 



C 8 -(3C 2 + C 1 )^( 



4m? m? . 



9 J 



mt 



+ 3ag 



LD 



(39) 



includes the contribution 



4, 4 8 8.., 1 , r - 

— my H — x H (1 H — xWx 

9 y 9 27 9 V 2 ; 



1 arctan 



=0(a;- 1) 



i(l + ix)v / T^(ln j^=§ + <7r)6l(l - x) 



(40) 
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from the operators 0\ 2 i and the long distance contribution 



9ld 



37T * , 

73 



a 



V=J/ip,ip 



m v T(V -> /+/-) 



(41) 



from the J / ip and ^' resonances. The parameter a that multiplies g^o is the 
same as in eq. || (for simplicity, I have taken the same parameter for both 
the J /ip and the tp' resonances, but this is not necessary); the relative sign 
between the long distance and short distance contributions is well determined 
[IHJ, up to the sign of a. The hadronic matrix elements are parameterized by 



< K\sYb\B > 
< K\sia^ (I + ^)b\B > 



{pB+PKYU{q 2 )+q»f-{q 2 ), 
s(q 2 )[(VB+VKYq v -(VB + VKfq 11 
+te^ap(pB+p K ) a q% (42) 



where, in the static 6-quark limit ||20|1 , s = — (/-+ 
gives 



/_)/2m j B- The modified 



BSW model 12 



with h Q 
then 



U(q 2 ) 

f-(q 2 ) 

0.379 and m + 



h 



1 



-/+(? 



Q+ (m B +m K ) 2 

2 m B - m K 



(43) 



5.89 GeV. The differential branching ratio is 



1 dT 

r dz 



^r B GW\v tb v; s \\^f{i- z ffl 

x(|C 9 | 2 + |C 8e// + 2C 7 | 2 ) 



(44) 



q 2 /m B ), where the lepton and kaon masses were neglected. This is 



shown in fig. 1 for a positive and negative. Studying the region of the 
interference between the short and the long distance contributions will allow 
to determine the sign of a, and resolve the ambiguity in eq. [T4|. At present, 
the necessary sensitivity has not been reached yet, and only an upper limit 
exists on the non-resonant B —± Kl + l~ decays ||2T . 
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3.2 "Non-factorizable" terms in B — > Dir 

The fortunate interference that allows to determine the sign of a is quite 
unique, and no similar effect appears for the decays of the type b — > cud, 
that would allow to determine the sign of a\ in eq. [15]. As for the final state 



interaction phase 5 in eq. |23|, it is known |22[ that it should be the same 
phase that appears in D-n elastic scattering, at the energy E cxn , — tub- But 
this is of little use in determining its value. Indeed, it is hard to think of 
an experimental test that would lift the 4-fold ambiguity in the values of the 
two "non-factorizable" terms in the B — > decays. On the other hand, it 
has been assumed throughout the analysis that the vacuum insertion result 
is a good approximation for the matrix elements of the operators with the 
correct color assignments. If it is further assumed that the vacuum insertion 
result should provide a first order approximation for the matrix elements of 
the color mismatched operators, then the A-terms should not be larger than 
unity. In particular, the solution for the B —> Dti decays is X\ ~ —0.17 and 



A 2 ~ 0.16, as in eq. 31 



Although there is no reason to expect that this is so (the arguments in 
ref. |J, for example, cannot be extended to the case of the "non-factorizable" 
terms), it should be pointed out that the solutions with small A-terms tend 
to agree with the values predicted by the theoretical calculations that are 
presently available. Using QCD sum rules techniques, it has been predicted 
that X\ ~ —0.33 (and, for the case of B — > J/4>K, X is between —0.30 
and —0.15 ||). Similar results have been obtained for the "non-factorizable" 
term in the amplitude for B° — B° mixing. The hadronic matrix element in 
the mixing amplitude is 

M mix . =< B° q \q a Y(l ~ 75)^7 M (1 - 75)M£° > • (45) 

Proceeding in a similar way as for the decay amplitudes, a "non-factorizable" 
term is defined by 

M mix . = 2(1 + — )m 2 B f B 

+ < B° q \q al »{l - ls)b a q^{l - ^)bp\W q > non . fact . (46) 

(where < B q {p)\q r y^{l — 75)6)0 >= —ifsp^ 1 )- The deviation from the vacuum 
insertion result (the first terms on the RHS of eq. 4"6) is parameterized by 
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the bag parameter Bb- Here, 

Bb — 1 + 7^ — j_, X mix ., (47) 



where 



< B° q \q a Y(l ~ l,)b a q^{l - 7 5 )M^° > 



"i/im \ - ioj <> i.) \- 10/ ■■II-— <j ' non-fact. fAQ\ 

mix. = q 2 f 2 ' \ / 



2 ^2 
B 



At present, the magnitude of _B# cannot be derived from the data on B® — B® 
mixing, because one lacks a precise determination of and of \V t d\ (for 
B° s — B® mixing, the CKM parameter is better known, but the strength 
of the mixing has not been determined experimentally). Using the lattice 
(B B = 1.2±0.2) H and the QCD sum rules (B B = 1.0±0.15) || estimates 



for the matrix element in the mixing, it follows that 

< X mix . < 0.53 and - 0.2 < X mix . < 0.2, (49) 

respectively. 



4 Conclusion 

The size of the "non-factorizable" terms in the amplitude for the decays 
B — > J/ipK and B —>■ Dn was derived from the experimental value of the 
corresponding branching ratios. The results can only be determined up to a 
discrete ambiguity. In the case of B — > J/tpK, the 2-fold ambiguity can be 
lifted by determining the sign of the interference between the short distance 
contribution to B — > Kl + l~, and the long distance contribution due to B — > 
KJ/ip — > Kl + l~ . A similar ambiguity appears in the case of the B — > 
decays; because the final state interaction elastic phase between the two 
isospin amplitudes can only be determined modulo ir, the ambiguity becomes 
4-fold. Contrary to the previous case, there is no simple way to determine 
the correct solution experimentally. However, all but one of the solutions 
indicates large "non-factorizable" terms that would indicate a breakdown of 
the vacuum insertion approximation when applied to the color mismatched 
operators. The analysis that was presented can be improved with future 
experimental results, in particular, with a measurement of the branching ratio 
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for — > D°7i°. Also, other 5-decays, similar to the ones shown in here, can 
be considered; the ambiguities in the values of the "non-factorizable" terms 
that are extracted from the data will appear in the same fashion. 
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Figure Caption 



Figure 1: Differential branching ratio for B — > Kl + l~~ (z = (p B — p K ) 2 /m 2 B ). 

The full line corresponds to the long distance contribution alone, 
whereas the other curves include the short distance contribution: 
with a > (dashed line) and a < (dotted line). 
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